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Abstract 

We consider some possible approaches to the fractional-order generalization of defini¬ 
tion of variation (functional) derivative. Some problems of formulation of a fractional- 
order variational derivative are discussed. To give a consistent definition of the fractional- 
order variations, we use a fractional generalization of Gateaux differential. 
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1 Introduction 

Derivatives of non-integer order are well-known in mathematics n-ia The fractional calculus 
has a long history from 1695, when the derivative of order a = 0.5 has been described by Leibniz 
[SI [71 [HI m [To]. Derivatives and integrals of fractional order have found many applications in 
recent studies in mechanics and physics. The interest to fractional equations has been growing 
continually during the last few years because of numerous applications El-ini- 

In mathematics and theoretical physics, the variational (functional) derivative is a gen¬ 
eralization of the usual derivative that arises in the calculus of variations. In a variational 
(functional) derivative, instead of differentiating a function with respect to a variable, one dif¬ 
ferentiates a functional with respect to a function. The fractional generalization of variational 
(functional) derivative was suggested in this paper. 

In this paper, some problems of formulation of a fractional-order variational derivative are 
discussed. To give a dehnition of fractional variation, we suggest to use a fractional generaliza¬ 
tion of Gateaux differential. 

In Section 2, some properties of the Riemann-Liouville or Gaputo fractional derivatives are 
noted. In Section 3, we give dehnitions of variational (functional) derivatives of integer orders. 
In Section 4, we discuss problems of different possible ways to dehne a fractional generalization 
of variational (functional) derivatives. In Section 5, we suggest a dehnition of fractional-order 
variational (functional) derivatives by using the proposed fractional generalization of Gateaux 
differential. In Section 6, a fractional variation of helds that is dehned by fractional exterior 
derivatives are considered. A conclusion is given in Section 7. 
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2 Fractional Derivative 


The theory of derivatives of non-integer order goes back to Leibniz, Liouville, Riemann, Griinwald, 
and Letnikov II1-0- The authors of many papers use the fractional derivative D'^ in the 
Riemann-Liouville or Caputo forms [1]. Let us give dehnitions of these derivatives and some 
properties. 


Definition 1. The Riemann-Liouville fractional derivative of the function f{x) belonging 
to the space AC"‘[a,b] of absolutely continuous functions is defined on [a,b] by the eguation 


D:f{x) = 


d” 


f{y) dy 


Vim - a) dx^ Ja (x - ’ 


( 1 ) 


where r( 2 ;) is the Gamma function, m is the first whole number greater than or egual to a. 


In equation ([T]), the initial point of the fractional derivative [T] is set to zero. The derivative 
of powers fc of x is 


Dfix - af = 




X 


(x > a), 


r(/c 1 — a) 

where k>l, and a > 0. Note that the derivative of a constant C need not be zero: 


( 2 ) 


X,- 




-c. 


a 


Therefore we see that constants C in the equation V (x) = C cannot dehne a stationary state 
for the equation DfV (x) = 0. In order to dehne stationary values, we should consider solutions 
of the equations Df,V{x) = 0. 

The Riemann-Liouville fractional derivative has some notable disadvantages in physical ap¬ 
plications such as the hyper-singular improper integral, where the order of singularity is higher 
than the dimension, and nonzero of the fractional derivative of constants, which would entail 
that dissipation does not vanish for a system in equilibrium. The desire to formulate initial value 
problems for physical systems leads us to the use of Caputo fractional derivatives [221 ESI ES] 
(see also mm rather than Riemann-Liouville fractional derivative. 


Definition 2. The Caputo fractional derivative of the function /(x) belonging to the space 
AC"'[a, 6] of absolutely continuous functions is defined on [a,b] by the eguation 


Dffix) = ^ r 

L [m — a) Ja 


dy d^f{y) 


(x > a). 


r(m - a) Ja (x - dy^ 

where f^^\y) = d^f{y)/dy'^, and m is the first whole number greater than or egual to a. 


( 3 ) 


This dehnition is of course more restrictive that ([T]), in that requires the absolute integra- 
bility of the derivative of order m. The Caputo fractional derivative hrst computes an ordinary 
derivative followed by a fractional integral to achieve the desire order of fractional derivative. 
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The Riemann-Liouville fractional derivative is computed in the reverse order. Integration by 
part of ([3]) leads us to 


ra—l ^k—a 

D:J(x) = D»J(x) - Y. Wj---fT/“’( 0 +). ( 4 ) 

^0 r(fc - a + 1) 

It is observed that the second term in Eq. (jH) regularizes the Caputo fractional derivative to 
avoid the potentially divergence from singular integration at a; = 0+. In addition, the Caputo 
fractional differentiation of a constant results in zero. If the Caputo fractional derivative is used 
instead of the Riemann-Liouville fractional derivative then the stationary values are the same 
as those for the usual case {y{x) — (7 = 0). The Caputo formulation of fractional calculus can 
be more applicable to definition of fractional variation than the Riemann-Liouville formulation. 


3 Variational Derivatives of Integer Order 

In mathematics and theoretical physics, the variational (functional) derivative is a generaliza¬ 
tion of the usual derivative that arises in the calculus of variations. In a variational (functional) 
derivative, instead of differentiating a function with respect to a variable, one differentiates a 
functional with respect to a function. 


3.1 Definition by increment and Taylor series 

The variational derivative can be dehned by the following way. Let us consider the increment 
of the functional F[u] that is dehned by the equation 

/\F[u] = F[u + h] - F[u], (5) 

and consider an integer-order variational derivative. 


Definition 3. If this increment of the functional F[u] exists, and ean he represented in the 
form 

AF[u] = 6F{u, h) + u{h,u), (6) 


where 


lim 

||/l||—40 


u{h,u)\\ 

\\h\\ 


0 , 


then 6F is called the first variation or Frechet derivative [FR ES of functional F. 
function h = h{x) is called the variation, and it is denoted by 6u. 


The 


Example 1. Let us dehne the functional 

rx2 

F[u] = / f{u)dx 
J Xi 


(7) 
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in some Bahan space E. The increment of the functional F[a] is defined by the equation 

/ X\ 


PX 2 rX2 

AF[u\ = F[u + h] — F[u] = / f{u + h)dx— f{u)dx. 

Jx\ Jx-[ 


( 8 ) 


Here we suppose that the variation h{x) = 6u{x) is equal to zero in boundary points xi and 
X 2 . Expanding the integrand f{x,u + h) in the power series up to hrst order 

d f {ij?\ 

f{u + h) = f{u) + ^ h + R{h, u). 


where 

Therefore we have 


du 


lim h ^R(h, u) = 0. 

h —^0 


sFiu]=r^hdx. 

J x\ UU 


The variational derivative for functional (JT]) has the form 

5F[u] df{u) 

6u du 

Example 2. Let us consider the functional 

rx2 

= / f{u,u^)dx. 

J Xi 

We can derive the hrst variation of the functional by the equation 

rX2 rX2 

5F[u] = F[u + 5u] — F[u] = / f{u + h^[uFh]x)dx — / f{u,Ux)dx. 

Jx^ J X^ 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Here we suppose that the variation h[x) = 5u{x) is equal to zero in boundary points Xi and 
X 2 - Let us expand the integrand f{u + 6u, {u + 5u)x) in the power series up to hrst order 

df df 

f{u + h,{u + h)x) = f{u, Ux) + -^h + -^hx, 

C/ ci C/ tix 

where hx = dh/dx. Then we get the variation of the functional 


rX2 

'df , df ' 

+ 7- — hx 

Jxi 

du dUx 


5F[u] = 

Integrating the second term by part, and supposing 

h{xi) = h{x2) = 0 , 

we get the result 

6F[u] = 


dx. 


(14) 


rX 2 

\df 

d 

(df\-\ 

J xi 

du 

dx 

\dux) _ 


h{x) dx. 


(15) 
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Using this relation, we get the variational derivative of the functional flT^ in the form 

SF[u] ^^_d_ 

5u du dx \duxj 

Example 3. If we consider the functional 

F[u] = f f{x,u,Ux,...u^'^^)dx 

Jxi 

then the variation of the functional is dehned by the equation 

^F[u] = J T.^iSuY^^dx, 

m=0 

and the variational derivative has the form 

dF[u] ^ ( df \ 


( 16 ) 


(17) 


(18) 


du ' dx^ J ’ 


(19) 


where = d'^u{x)/dx^. 

3.2 Definition by composite function of parameter 

There is the other approach to dehnition of variation and variational derivative. Let us consider 
the functional flT^ . where m is a function of coordinates x and parameter a, i.e., u = u{x,a), 
and 


F\u] = 


r ^2 


f{u{x, a),Ux{x, a))dx. 


( 20 ) 




The derivative of F with respect a can be derived 


dF r*2 ^ df r ^'2 
dd J XI dd J x\ 


df du ^ df dux 


du dd dux dd 


dx. 


( 21 ) 


Using 


the conditions Su{xi,d) 


dux d du d du 
dd dd dx dx dd ’ 

6 u{x 2 , a) = 0, and integrating by part, we get 


As the result, we have 


dF 

/■*2 du 

'df 

A 


dd 

Jxi dci 

du 

dx 

\dux) _ 


dF du{x) 6 F[u]^ 

dd Jxi dd 6 u{x) 


( 22 ) 


(23) 
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A fractional-order generalization of this approach is difficult to realize. The reason of this 
difficulty is connected with difficulty to generalized the rule of differentiating a composite 
functions 

dfiuia)) df du 

V = 24 

da ou oa 

for the fractional case. It is known that coordinate transformations is connected with the 
derivative of a composite function Dlf{u{x)) = {Dl^f){u = u{x)) {D].u){x). The formula of 
fractional derivative of a composite function (see Equation 2.209, Section 2.7.3 page 98 i) is 
the following 


k=i a -I- ij 


1 /(T)»(x))xa, 


r(l -a)' 


r=l ^r- 


T\ 


(25) 

where the sum X) extends over all combinations of non-negative integer values of Oi, 02 , . . . , 
Ofc such that 


rar = k, Y^r = 


m. 


(26) 


r=l 


4 Problems to Formulate Fractional-Order Variational 
Derivative 

Let us consider some problems in the formulation of the fractional-order generalization of 
variational derivatives. 

Problem 1. To dehne the fractional variation, we can use the fractional Taylor series (see 
Section 2.6 in m)- Expanding the integrand f{u + h) in the fractional Taylor power series 

/(«■ + 1 ‘) = \ (hr D:f(u) + _-L_(ft)”+‘£>;/(t.) + .... (27) 

1 (a -I-1) 1 (a -I- 2) 

we can see that this series cannot have the term f{u). Therefore, we cannot consider the 
increment 

^fiu) = f{u + h)-f{u). (28) 

In order to use Taylor series for the dehnition, we can consider the increment in the form 

= f{u +h) - f{u +qh). (29) 

where 0 < g < 1. In this case, we have 

= n + n + .... (30) 

r(a-I-1) r(Q;-|-2) 

This approach can be connected with the q-analysis [3l], and fractional q-derivatives [321133] . 
We will consider the fractional variational q-derivative in the next paper. 
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We can consider the increment of the fnnctional F[u] that is defined by the equation 

AqF[u] = F[u + h]— F[u + qh]. (31) 

If this increment of the functional F[u] exists, and can be represented in the form 

AqF[u] = S^F{u, h") + u{h^, u), (32) 


where 


lim 

|| fe||—>-0 


|a;(h", u) 

~w\r 


= 0 , 


then S°‘F can be considered as a fractional variation of functional F. 


Problem 2. To define the fractional generalization of variation and fractional exterior 
variational calculus [21], we can use an analogy with definition of fractional exterior derivative. 
If the partial derivatives in the definition of the exterior derivative 


d = dxid/dxi 


are allowed to assume fractional order, a fractional exterior derivative can be defined [20| by 
the equation 

o!“ = (dxXDl, (33) 

where are the fractional derivative with respect to x. Using this analogy, we can dehne the 
fractional variation by the following way. For the point u of functional space, we can define the 
fractional variation 6 F[u] of the functional 

rx 2 

= / f{u,u^)dx, 

J Xi 


where Ux = du/dx, by the equation 


6 ^ Flu] = 




dx 


{Su)“D“f(u,u^) + (6u,)“D“J(u,u, 


(34) 


Unfortunately, this approach leads to the difficulties with the realization of integration by part 
in the second term. It is easy to see that the variation {6ux)‘^ cannot be represented as some 
operator A{d/dx) acts on the variation (hu)", i.e., we have 


d 


d 


{duxix))"" = —6u{x) ^ Al — ] (5n)“. 


dx 


dx, 


(35) 


In the particular case, {6ux{x))°‘ = {d{6u)/dx)°‘ ^ d{6u)'^/dx if a 7^ 1. 

Problem 3. To define the fractional variational derivative, we can use the fractional 
derivative with respect to function (see Section 18.2 in [I])- The fractional Riemann-Liouville 
dehvative of the function f{x) with respect the function u{x) of order a (0 < a < 1): 


Do+,uix)fix) 


1 / du{x) \ ^ r ^ f{y) du{y) 

r(l —a) \ dx ) dx Jo ^ [u{x) — u{y)]'^ dy 


(36) 
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Using the function u = u{x,a), we can define fractional derivative with respect to function of 
parameter 

-1 


na ft \ ^ du{x,a) 


- 7 - / dbr 


f{x,b) 


du{x, b) 


da J da Jo [u{x, a) — u(x, b)Y 
As the result, the fractional variation can be dehned by the equation 


db 


rx2 

6 ^F[u{x)]= / dxDY(^^^^-^f{u{x,a)){6u{x,q))'^, 

f X\ 


5^F[u] 


= -D“(i,a)/(«(a:))' 


(37) 


(38) 


This approach leads to the difficulties with the realization of integration by part in the second 
term fl35l) . In this paper, this approach to definition is not considered. 

To avoid all these difficulties and problems, we suggest to dehne a fractional variational 
derivative by using a fractional-order generalization of Gateaux differential. 

5 Fractional Generalization of Gateaux Differential 

In this section, we consider 3 steps to define the fractional generalization of variation by using 
some generalization of Gateaux differential. We consider these steps in order to explain the 
final dehnition of fractional Gateaux variation. 


5.1 Variations of integer order 

Suppose the functional F[u] is continuous (smooth) map (with certain boundary conditions) 
from everywhere dense subset D{F) of Banah space to space R. Let us define the Gateaux 
differential [28l |29l |30] of a functional F[u] at the point u{x) of subset D{F) of the functional 
Banah space. 

Definition 4. The Gateaux differential (or first variation) of a fucntional F[u] is defined 
by the equation 


5F[u, h] = i^F[u + eh] j = lirn 


F[u -I- e h] — F[u] 


(39) 


e=0 


if the limit exists for all h{x) G D{F). The function h{x) is called a variation of function u{x) 
and denoted by Su{x) = h{x). 

A first variation of functional F[u] at the point u = u(x) is defined as a first derivative of 
functional F[u + eh] with respect to parameter e for e = 0. 

If the Gateaux’s differential of the functional 


F[u] = J dxf{x,u,u ^,...) 
is a linear with respect to h{x), then we can write 

SFlu,h] = jda:E(.,u,u^,...Ma:), 


(40) 


( 41 ) 







and E{x,u,Ux, ■■■) is called the variational derivative and is denoted by 6F/6u. 

Example. For the example, we can consider the functional 

F[u] = / [u{x)Ydx. (42) 

J Xi 

The functional F[u + eh] has the form 

rX2 rX2 rX2 

F[u + eh] = / [u{x)]‘^dx + 2e u{x)h{x)dx + / [h{x)f‘dx. (43) 

J Xl J XI J XI 

The derivative of first order with respect to parameter e is 

d /'^2 PX2 

—F[u + eh]=2 u{x)h{x)dx + 2e [h{x)]‘^dx. (44) 

(!€ J x\ J x\ 

As the result, the variational derivative is equal 

C d \ 

—F[m + eh]\ = 2 / u{x)h{x)dx. (45) 

de J Jxi 


5.2 STEP 1. 

It seems that we can define a fractional variation by the equation 


6 ^F[u,h] = 


d 

de 


F[u + eh] 


(46) 


e=0 


Let us consider the fractional Riemann-Liouville derivative of functional fl43|) with respect 
to parameter e: 




fX2 rX2 PX2 

F[u + eh] = (Dfl) / [u{x)]‘^dx + 2{D°e) / u{x)h{x)dx + {D^e'^) / [h{x)]‘^dx. (47) 

Jxi J x^ J x^ 


Using equation 


T^a k ^ r(fc + l) 


r(fc + 1 — a) 


(48) 


we have 


+ eh] = 


[u{x)f‘dx + 


2 e 


1 — OL 


r(l — a) Jxi ^ r(2 — a) Jxi 

9 rx 2 

' / [h{x)]^dx. 


u{x)h{x)dx+ 


r(3 — a) 


(49) 


fX\ 


If the derivatives with respect to parameter e in the definition of variational (l39|) are allowed 
to assume fractional order, a fractional variational can be defined. Unfortunately, if we define 
the fractional variation of the functional by the equation 




(50) 
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then we have some problems that state about incorrectness of definition fl50l) . These problems 
are following: 

1) The hrst term of right hand side of equation fl49l) leads us to infinity. If e tends to zero, 
then we get —)■ oo, and D'^F[u + eh] —)■ cx3. Therefore the hrst term that is follows from the 
relation D'^C = 0 must be removed. For this aim we can use the Caputo fractional derivative 
[21E21 ESI El]- Using (H]) for 0 < a < 1, we have 

D- J{e) = D:fie) - ^,^1^/(0+). (51) 

2) The second term of the right hand side of equation fH9|l is proportional to This 

proportionality leads us to zero in the limit e —)■ 0, and we cannot derive some nonzero relation. 

Therefore we must consider the functional F[u + (eh)“] in the dehnition. 

As the result, we cannot use the dehnition fll6l) and fl50H . 

5.3 STEP 2. 

Let us consider the following equations for the dehnition of fractional variational functional 

F[«]: 

rFH = (£>fF|«+(efc)"])„„. (52) 

where is a Caputo fractional derivative of order a with respect to e. 

If we consider the functional flT2]) . then 

rX2 rx2 rX2 

F[u + {ehY]= / [u{x)Y‘dx + 2e'^ j u{x)h°'{x)dx + e‘^‘^ [h{x)]‘^°'dx. (53) 

Jxi jXI jXI 

The fractional derivative fl5^ of this functional is equal to 

DfFlu + (eh)"] = 2r(a + 1) / u{x)h°‘{x)dx + f \h[x)]‘^^dx. (54) 

Jxi r(Q; + l) Jxi 

Therefore, we get the following equation 

rx2 

(I1"F[m + (eh)"])^^Q = 2r(a + 1) / u{x)h°'{x)dx, (55) 

J Xi 

Note that for a = 1, we get the usual relation 

rX2 

6°‘~^F[u] = 2 / u{x)h{x)dx. (56) 

Jxi 

Unfortunately, if we use a = 0, we get 

rx2 rx2 

F[u\ = = 2 / u{x){h{x))^dx = 2 / u{x)dx ^ F[u]. (57) 

Jxi J X\ 

Therefore we cannot use the dehnition ([52]). 
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5.4 STEP 3. 

The usual definition 


can be rewritten in the form 


5F[u] = 


de 


n 1 + 


eh' 


u 


(58) 


e=0 


The fractional variation can be defined by the equation 

eh 
u 


S°F[u] = D'^F 


n 1 + 


(59) 


e=0 


This definition is more consistent in order to realize the physical dimensions. In this definition 
the parameter e is dimensionless. 

If we consider the functional flT2|) . then 

^ =F[u+ (eh)“n^-“] = 

rX 2 rX 2 rX 2 

= [u{x)]‘^dx + 2e‘^ u^-'^{x)h'^{x)dx + e^‘^ u^-^^{x)[h{x)Y'^dx. (60) 

J XI Jxi Jxi 

The fractional derivative of this functional is equal to 



eh 


1 - 

s 

+ 

u 

)J 


PX2 

D’^F[u + (eh)“n^-“] = 2T{a + 1) / u^-'^{x)h‘^{x)dx + 

J Xl 


r(2a + l)^„ r\.2- 


r((T -|-1) Jxi 

As the result, we get 

(D“F[n + = 2r(a + 1) u^-^{x)h^{x)dx 


rX 2 

/ u^-^'^{x)[h{x)f^dx. 

J X^ 

(61) 


For a = 1, we get the usual relation. For a = 0, we have 

fX 2 
/ X\ 

If we consider the functional 


PX2 

F[u] = (5“=°F[m] =2 u^{x)dx ^ F[u]. 

J X-\ 


rx 2 


F[u] = I u'^{x)dx, 

J Xl 

then definition fl5^ leads us to the relation 

(^DfF[u + (eh)"n^-"])^_^ = nT{a + 1) u^-^{x)h^{x)dx 


and, for a = 0 


PX2 

F[u] = (5“=°F[n] = n / u^{x)dx ^ F[u]. 
J X^ 


(62) 

(63) 

(64) 

(65) 

( 66 ) 


Therefore we must realize some changes in the definition (159]). This modification is suggested 
in the next subsection. 
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5.5 Definition of variation of fractional order 

Taking into account the remarks in the form of Step 1-3 and using Eq. fl58|l for variation of 
integer order, we can define the fractional variation in the following form. 


Definition 5. The fractional-order variation 5°'F[u] of the functional F[u] is defined by the 
equation 

= (^DfF 

or, in an equivalent form 



M 1 + 


eh' 


u 


J / e=0 


(67) 


S‘‘F\u] = [d“F + 


'e=0 


( 68 ) 


where e > 0. 


Example. If we consider the functional fH2|) . then we have 


F 





u^-^'^h^^^ie + u/hf^dx. 


(69) 


Using the relation 

Dnc + u/h)^ = t^^±f-(e+u/hr, (70) 

i (a 1) 

where we consider e -|- u{x)/h{x) = e — cq with cq = eo(x) = —u{x)/h{x), we get the fractional 
derivative of functional (1691) in the form 


D<fF 


u 


1 — 0. i 


u + ehy 


r( 2 a + 1 ) 

r(Q( -|-1) 





-|- u/h)^dx. 


As the result, we get 





e=0 


r( 2 a + 1 ) 

r(o( -|-1) 



dxu^-‘^{x)h‘^{x). 


(71) 


(72) 


For 0 = 1, and o = 0, we get the usual relations. Therefore definition fl57)) satishes the cor¬ 
respondent requirement for o = 0 and o = 1 As the result, we get the fractional variational 
derivative of order 0 < a < 1 in the form fl68|) . 


Proposition 1. The fractional-order variation derivative ( [h7| ) of the functional 

rX 2 

F[u] = / u'^{u)dx. 

Jxi 


has the form 

rX2 

5°F[u] = X{a,n) / dx {Dy^){x)h^{x), 

Jxi 


(73) 

(74) 
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where 


A(a, n) 


T{na + l)r(n + 1 — a) 
r((n — l)a + l)r(n + 1) ’ 


(75) 


Proof: Using fl73|) . then we have 


F 


+ ehy 


rX2 


u^-^^h'^^ie + u/hY^dx. 


IXl 


Using the relation 




(76) 


(77) 


we get the fractional derivative of functional fl75]) in the form 

6 ^F[u] = (d^F [u^-^{u + dxu'^-^{x)h^{x). (78) 

This equation can be represented in the form 

rx 2 

8 ^F= / dx{Dy^)h^{x), (79) 

J XI 

up to numerical factor. For equation (1791) . we have 

( 0 >") = (80) 

1 (n + 1 — a) 

This end of proof. 


The fractional variational derivative of order a > 1 can be dehned by the usual relation 

5“ = (81) 


by analogy with fractional derivative 


D 


a 

X 


do:!"] 


Di'^K 


(82) 


where [a] is a whole part of a, and {a} is a fractional part of number a, i.e., {a} = a — [a]. 


5.6 Functional with derivative of field 

Let us consider the functional F[u] that has the form 

fX2 

-^N = / f{u,U:^)dx (83) 

J Xi 
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where Ux = ■ The hrst variation of this functional is defined by 

ax 




e=0 


C d /■^2 \ 

— j f{u +eh,[u +eh]x) j dx, 


where we use 


In order to use the equation 


[u + e/z]i ^ 


(84) 

(85) 


Bf(e - to)'’ = - fo)"-". (86) 

r(/f + 1 - a) 

for the relation [e + u{x)/h{x)], we consider eo = eo((c) = —u{x)/h{x). The dehnition of the 
fractional variation of functional fl53D can be realized in the following form. 


Definition 6. The fractional-order variation derivative of the functional l[83) is defined 
by the equation 

S°F[u] = lim (D“/+ eh)“, [u^-^{u + eh)%)) . (87) 


Using the definition fl 68 D . we have the fractional variation of functional fl83l) in the form 

S^F[u] = [Dff(u^-^{u + eh)", [u^-^{u + . ( 88 ) 

It is easy to see that the expression [u^~°‘{u + eh)°‘]x cannot be represented in the form of similar 
Eq. (|85il . This expression has the form 

[m^“"(m + eh)^]x = [(1 — a)u~^{u + eh)" + au^~°‘{u + eh)"“^]Ma, + au^~°‘{u + eh)°‘~^ehx. (89) 
For a = 0, we have 

[m^""(m + eh)"]^, = Ux, (90) 

and for a = 1 , we have 

[u^~°‘{u + eh)"]a, =Ux + ehx- (91) 

Let us give the proposition for a special form of the functional. 

Proposition 2. The fractional-order variation derivative |g7| ) of the functional 

rx2 

F[u] = / u{x)ux{x)dx. (92) 

J Xl 

has the form 

S°‘F[u]= ( dx (^Ai{a)ux{x)h{x) + A 2 {a)u^~'^{x)ah^~^{x)hx{x)^ , (93) 
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where 


AAo) = 


(1 — a)T(2a + 1) aT{2a) 


r(Qf +1) 

Proof: Using the functional 


r(«) 


Aoioi) — 


o^r( 2 o) 
r(o + 1 ) 


PX2 

+ eh)"] = / dx = + eh)°'[u^~°‘{u + eh)"]^; = 

J X\ 

rx2 / 

Jxi ^ 

+ aM2-2a;^2a-i ^ u/hf^-^eh^). 

Let us consider cq = eo(x) = —u{x)/h{x). We can use the following relations 

“r(/3+l-a)’ 


(94) 


(95) 


(96) 


and 

Df (6 - 6 o)'’(e - 0)-- = (-7. /3 + 1; )J + 1 - a; -^) («o)’(£ - (97) 

Note that Eq. is satisfyed for [5 > —1, and Eq. is satisfyed for [5 > —1, e > Cq > 0 . 
Here F 2 -i{a,b,c, z) is the Gauss hypergeonietric function [T]: 


F2;i(a,6, c,2; 


OO 


) = E 

k=0 


{a)k{bk) 

{c)k kV 


and 


{z)k = z{z + 1 )...( 2 : + n - 1 ) = 


r(z + n) 


For the function F 2 ;i(a, b, c, z) exists the Euler representation 

r(c) 


F2;i(a,6,c,z) = —^(1-^)' '' ^(1-^^) ‘"dt, 


V{b)V{c-b) 

where Re{c) > Re{b) > 0, |ar 5 f(l — z)\ < n . Using 

r r j, 11 r(c)r(c-a- 6 ) 

^ 2;1 a, b; c; 1 = —-—-- 

r(c — a)r(c — b) 

where Re{c — a — b) > 0 . For the case a > 1, we can use 


(98) 


(99) 


17 r 1 («)fc(^fc) r r , 7 , 7 ,71 

■^F 2 -i[a, b] c; z\ = —- F 2 -i[a + fc, 6 + fc; c + fc; ^J, 


( 100 ) 
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where k = [a] is a whole part of a. As the result, we get 


r ^2 


6 ^F[u] = / dx 

J Xl 


+ 


(1 —a)r(2a + l) ar(2a) 
r(a + l) ^ r(a) 

'x\ r(tt +1) 


Ux{x)h{x)+ 


rx 2 


(101) 


This end of the proof. 

For a —)• 1, we have the usual relation 

6 F[u] = I dx[uxh + uhx\. 

Note that we can use 


fX 2 

> Xl 


= {h°‘)x, 


and then 




da:—T—Tvrt- tU^ °‘ah°' ^hx = 




( 102 ) 


Jxi r(—a)r(a; + 1) Jxi r(—Q;)r(a + 1) 

This allows us to realize the integration by part in second term of Eq. fllOip . 

6 Fractional Variation of Fields 

In this section we consider the dehnition of fractional variational derivative without using the 
increment 

AF[u] = F[u + h]- F[u] 

of the functional F[u], and without using the derivative with respect to parameter e as in the 
Gateaux derivative. We suppose that functional F[u\ has some densities f{u,Ux ,...). 

If u = (ui, be a smooth vector-function that is defined in the region W C i?", 

then the variation of the functional 

F[u] = / f{u,Ux)dx 

Jw 

can be defined by the relation 

6 F[u]= [ 6f{u,Ux)dx= f 
Jw Jv 


5u^ + -rArSut 




dux 


dx. 


(103) 


To define the fractional generalization of variation and fractional exterior variational calculus 
izil. we can use an analogy with dehnition of fractional exterior derivative. If the partial 
derivatives in the dehnition of the exterior derivative 


d = dxid/dxi 
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are allowed to assume fractional order, a fractional exterior derivative can be defined 1211 by 
the equation 

r = (104) 

where D" are the fractional derivative with respect to x. Using this analogy, we can define the 
fractional variation by the following way. For the point u of functional space, we can define the 
fractional variation 6F[u] of the functional 

r^2 

-^N = / f{u,u^)dx, (105) 

J X\ 


where = du/dx, by the equation 


5°‘F[u\= / 6f{u,Ux)dx= / dx 

JXi JX2 


{6u)^D^f{u,Ux) + {6ux)'"D^J{u,Ux) 


(106) 


This approach has a difficulty with the realization of integration by part in the second term 
of (llOOp . It is easy to see that the variation (hwa;)" cannot be represented as some operator acts 
on the variation (5 m)", i.e., we have 


In order to solve this difficulty we can use the following. 
Let us define the fractional variation of the functional 


by the equation 


by analogy with 


In this case, we have 


It is known, that 


F[u]= J f{Ui,U2,...,Um)dx, 


5"F[m] = J 


Y.iD:j)i6u,y 

.k=l 


dx 


= / f2iDlj)6ukdx. 

k=i 




Y,(DZ,n)(Su,Y 


.k=l 


Su, = j 




.k=l 


dx = J [Ski5{y - x)5uk] dx. 


Using the Caputo fractional derivative, we have 


Di,(y)Ui(x) = - x). 


(107) 

(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 

(113) 
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Here the Caputo fractional derivative leads us to the 6ki, i.e., 


Du,,{y)Mx) = 0 k^l. 


As the result, we get 


Therefore, we have 


S^ui = 


u 


1 — OL 


r( 2 -a) 


5ki{8ukT. 


(114) 

(115) 


(^Mfc)" = r (2 — a)u^ 

For the Riemann-Liouville fractional derivatives, we can derive some analogous relation: 

{5ukT = Am (r (2 - a)<-^ 1(1 - a)u^\r) ( 116 ) 

Substituting Eq. flll5p in Eq. (11091) . we have the following dehnition. 

Definition 7. The fractional order variation d°'F[u] of the functional MO^) is defined by 
the equation 


S“F[n] = r(2 -a) j 




.k=l 


dx. 


(117) 


We can dehne the fractional variation for derivative duk{x)/dx by the equation 

S^^Ukix) = ■^6°‘uk{x), (118) 

dx dx 

where we suppose that the helds Uk do not connected by some constraint. Analogously, we 
have 

6^D^Uk{x) = D^J^Uk{x). (119) 

Let us consider ui = m, and U 2 = Ux and the functional 

PVA = J f(,Ui,U 2 )dx = J f{u,Ux)dx. (120) 

As the result, we have 


6^F[u] = r(2 -a) J 




u 




f)u 


a—1 

kx 


5°‘ukdx. 


( 121 ) 


We can use this relation only if Ux and u can be considered as independent values, i.e., for the 
case 

DZu = t), DX = 0. (122) 

Note that Eq. (11221) cannot be satished in general case. Therefore, we must consider the 
functional 

F[t,] = j 


f{Ui,U2) + A f «2 - -^Ui 


dx. 


(123) 
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instead of the functional (11201) . 


Proposition 3. The fractional-order variation equation 

6^F[u] = 0 

of the functional M23i) gives 

(D."/)«?"-" - ^ hr'-r'fl"/) = 0 . 

where U 2 = dui/dx. 

Proof: The fractional-order variation of the functional (I123p gives 


6^F[u] = r(2 -a) j dx 


d 


P“/K- 


OL—l 




-Fr (2 -a) J dx [{D^J)u2 ^ + X(D^^U2)u^ ^"^2+ 

-I- r(2 — a) J dx 


U2 - -j-ui 
dx 


{DfX)Xf-^6^X. 


Using 


m) - (a<"^) , 


dx^ ^ U “1 ^dx 

we get the following fractional variation of the functional 


i»F|u] = r(2 -a) fdx (DZJ)ur' + 


5°‘uiF 


+r(2-a) dx {DU)ur^ + X{D:u 2 )ur^ S^U 2 + 


+ r(2 — a) j dx 

As the result, we get the held equations 


d 

U2 - -r^i 
dx 


(B“/)«r‘ + B”mT(A<-‘) =0, 

(D“/) + A(D“ «2) = 0, 

d 

U2 - r-Ul - 0 

ax 


(r>JA)ArM“A. 
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(124) 


(125) 


(126) 


(127) 


(128) 

(129) 

(130) 




















From the Eq. fll29p ,we derive the Lagrange multiplier 



(131) 


Substituting this equation in Eq. fll28p . we have 



0 . 


(132) 


Using D^u = "/r(2 — a), we get (11241) . 

Equation (11241) is the fractional Euler-Lagrange equation. 

7 Conclusion 

In the general case, the equation of motion cannot be derived from the stationary action 
principle. The class of equations that can be derived from stationary action principle by using 
fractional variation is a wider class than the usual class equations that can be derived by usual 
(integer, first order) variation. The usual equations of motion can be considered as special case 
of equations that can be derived by fractional variation such that a = 1. 

It seems that the fractional variations are abstract and formal constructions. For this reason, 
we would like to pay attention that suggested fractional variations can have the wide application 
in study of fractional gradient type equations and fractional generalization of Lyapunov direct 
(second) method in the theory of stability. 

The possible importance of fractional variations are connected with following ideas. The 
class of gradient dynamical systems is a restricted class of all dynamical systems. However these 
systems have the important property. The gradient system can be described by one function - 
potential, and the study of the system can be reduced to research of potential. For example, the 
way of chemical reactions is defined from the analysis of potential energy surfaces [aslIMlET]. 
The fractional gradient systems has been suggested in Refs. The fractional gradient 

systems are non-gradient dynamical systems that can be described by one function - some 
potential. For example, the Lorenz equations and Rossler equations are fractional gradient 
systems [HI [19]. Therefore the study the some non-gradient system can be reduced to research 
of potential. For example, the way of some chemical reactions with dissipation, dynamical 
chaos and self-organizing can be considered by the analysis of some potential energy surfaces. 
The suggested fractional variations allows us to define the fractional generalization of gradient 
type equations that can have the wide application for the description dissipative structures 
[381139] . The suggested approach can also be generalized for lattice systems by using the lattice 
fractional calculus [34] . 


20 












Appendix 

There are the following rules for variational dehvatives. 


The variation derivative of the field u{x) is defined by the equation 

5u{x) 


5u{y) 


= 5{y-x), 


where we use 


5u{x) = J 6{y — x)6u{y). 
The variational derivatives of linear functional 

F[u] = J g{x)u{x)dx 
can be calculated by the simple formula 


6F[u] 

Su{y) 


= j = j !l(x)S(y - x)dx = h(y). 


Su(y) 


(133) 

(134) 


(135) 


• If the u function in the functional is affected by differential operators, then, in order to 
make use of the rule (I133p . one should at hrst ’’throw them over” to the left, fulhlling 
integration by parts. For example, 

J g{x)[Vu{x)]dx = J[^g{x)]u{x)dx = -Vg{y). 

We assumed here that on the boundary of integration domain the product u{x)h{x) 
becomes zero. 


The variational derivative of nonlinear functionals is calculated according to the rule of 
differentiating a complex function similarly to partial derivatives: 




6u{y) 
For example. 


f{u{x))dx = 


6f{u) 5u{x) 
6u{x) Su{y) 


dx = 


-d[y — xjdx = 


6u{x) 


6u{y) 


6u{y) 


J[u{x)]"'dx = n[u{y)]"' ^ 
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